(N 
O 
O 
(N 

-t— > 
o 

o 

o 

m 



I 

O 

o 



> 
oo 

m 

o 

O 



I 

O 

o 



X 



Current- voltage characteristics of molecular 
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Abstract — This paper addresses the question of whether a 
"rigid molecule" (one which does not deform in an external 
field) used as the conducting channel in a standard three- 
terminal MOSFET configuration can offer any performance 
advantage relative to a standard silicon MOSFET. A self- 
consistent solution of coupled quantum transport and Pois- 
son's equations shows that even for extremely small chan- 
nel lengths (about 1 nm), a "well-tempered" molecular FET 
demands much the same electrostatic considerations as a 
"well-tempered" conventional MOSFET. In other words, we 
show that just as in a conventional MOSFET, the gate oxide 
thickness needs to be much smaller than the channel length 
(length of the molecule) for the gate control to be effective. 
Furthermore, we show that a rigid molecule with metallic 
source and drain contacts has a temperature independent 
subthreshold slope much larger than 60 mV / decade, because 
the metal-induced gap states in the channel prevent it from 
turning off abruptly. However, this disadvantage can be 
overcome by using semiconductor contacts because of their 
, band-limited nature. 

' Keywords — Molecular electronics, MOSFETs, elec- 
trostatic analysis, quantum transport. Non-equilibrium 
Green's function (NEGF) formalism. 



I. Introduction 

MOLECULES are promising candidates as future elec- 
tronic devices because of tlieir small size, chemical 
tunability and self-assembly features. Several experimen- 
tal molecular devices have recently been demonstrated (for 
a review of the experimental work see Q). These include 
two terminal devices where the conductance of a molecule 
coupled to two contacts shows interesting features such as 
• a conductance gap , asymmetry [|j and switching Q . 
[ Molecular devices where a third terminal produces a nega- 
' five differential resistance |^ , or suppresses the two termi- 
\ nal current |^ have been theoretically studied, but most of 
' the work on modeling the current-voltage (IV) character- 
' istics of molecular conductors has focused on two-terminal 
"devices (see, for example, §, g, [§, Q, (n), (l| and 
\ references therein) . 

The purpose of this paper is to analyze a three-terminal 
molecular device assuming that the molecule behaves es- 
sentially like a rigid solid. Unlike solids, molecules are ca- 
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pable of deforming in an external field and it may be pos- 
sible to take advantage of such conformational effects to 
design transistors with superior characteristics. However, 
in this paper we do not consider this possibility and simply 
address the question of whether a "rigid molecule" used as 
the conducting channel in a standard three-terminal MOS- 
FET configuration can offer any performance advantage 
relative to a standard silicon MOSFET. 

Although rigorous ab initio models are available in the 
literature ||] , , 1 11 1 , Q , they normally do not account 
for the three-terminal electrostatics that is central to the 
operation of transistors. For this reason we have used a 
simple model Hamiltonian whose parameters have been 
calibrated by comparing with ab initio models. We believe 
that a simple model Hamiltonian with rigorous electrostat- 
ics is preferable to an ab initio Hamiltonain with simplified 
electrostatics since the essential physics of a rigid molecular 
FET lies in its electrostatics. 

The role of electrostatic considerations in the design of 
conventional silicon MOSFETs (with channel lengths rang- 
ing from 10 nm and above) is well understood. For the gate 
to have good control of the channel conductivity, the gate 
insulator thickness has to be much smaller than the channel 
length. Also, for a given channel length and gate insulator 
thickness, a double gated structure is superior to a single 
gated one, simply by virtue of having two gates as opposed 
to one. If a molecule is used as the channel in a stan- 
dard three-terminal MOSFET configuration, the effective 
channel length is very small - about 1 nm. Would similar 
electrostatic considerations apply for such small channel 
lengths? In this paper we answer this question in the affir- 
mative. Specifically we will show that: 

• The only advantage gained by using a molecular conduc- 
tor for an FET channel is due to the reduced dielectric 
constant of the molecular environment. To get good gate 
control with a single gate the gate oxide thickness needs to 
be less than 10% of the channel (molecule) length, whereas 
in conventional MOSFETs the gate oxide thickness needs 
to be less than 3% of the channel length . With a dou- 
ble gated structure, the respective percentages are 60% and 
20% 11. 

• Relatively poor subthreshold characteristics (a temper- 
ature independent subthreshold slope much larger than 
60 mV/ decade) are obtained even with good gate control, 
if metallic contacts (like gold) are used, because the metal- 
induced gap states in the channel preclude it from turning 
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Fig. 1 

(a) Schematic of a Phenyl Dithiol molecule coupled to gold 

SOURCE and drain CONTACTS. A THIRD (GATE) TERMINAL 
MODULATES CONDUCTANCE OF THE MOLECULE. ThE PHENYL RING IS 
SHOWN IN THE PLANE OF THE PAPER FOR CLARITY, IN THE ACTUAL 
SIMULATION THE PHENYL RING IS PERPENDICULAR TO THE PLANE OF 
THE PAPER, "facing" THE GATE ELECTRODE, (b) ThE MOLECULE IS 
DESCRIBED BY A HaMILTONIAN H AND A SELF-CONSISTENT 
POTENTIAL Use- THE EFFECT OF THE LARGE CONTACTS IS 
DESCRIBED USING SELF-ENERGY MATRICES Si, 2- SCATTERING 
PROCESSES MAY BE DESCRIBED USING ANOTHER SELF-ENERGY MATRIX 

Sp. The source and drain contacts are identified by their 
RESPECTIVE Fermi levels /xi and ij,^- Given H, Use, ^i,2,p and 
Ml, 2 THE Non-Equilibrium Green's Function (NEGF) 

formalism has clear PRESCRIPTIONS TO OBTAIN THE DENSITY 
MATRIX FROM WHICH THE ELECTRON DENSITY AND CURRENT MAY BE 
CALCULATED. At EQUILIBRIUM (ZERO DRAIN BIAS) = fl2 = Ef, 
WHERE Ef IS THE COMMON EQUILIBRIUM FeRMI LEVEL OF THE 
CONTACT-MOLECULE-CONTACT SYSTEM. WHEN A DRAIN BIAS VuS 'S 

APPLIED THE SOURCE AND DRAIN FeRMI LEVELS SEPARATE BY AN 
AMOUNT EQUAL TO qVoS {Q- ELECTRONIC CHARGE) AND A CURRENT 

FLOWS. This non-equilibrium situation may be modeled by 

SELF-CONSISTENTLY SOLVING THE COUPLED NEGF-POISSON'S 

equations as shown in (c). The word "Poisson's" is in quotes 

AS A reminder that MORE SOPHISTICATED THEORIES LIKE THE 

Hartree-Fock or the Density Functional theory may be used 
TO obtain Use- 



off abruptly. Preliminary results with a molecule coupled 
to doped silicon source and drain contacts, however, show a 
temperature dependent subthreshold slope (~ fc^T/g). We 
believe this is due to the band-limited nature of the silicon 
contacts, and we are currently investigating this effect. 

Overall this study suggests that superior saturation and 
subthreshold characteristics in a molecular FET can only 
arise from novel physics beyond that included in our model. 
Further work on molecular transistors should try to capi- 
talize on the additional degrees of freedom afforded by the 
"soft" nature of molecular conductors - a feature that 
is not included in this study. 

Although there has been no experimental report of a 



moleculer FET to date^ , judging from the historical devel- 
opment of the conventional silicon MOSFET, it is reason- 
able to expect that a single gated structure would be easier 
to fabricate than a double gated one. With this in mind, 
in this paper we mainly focus on a single gated molecular 
FET geometry (see Fig. p. Few key results with a double 
gated geometry will be shown wherever appropriate to em- 
phasize the differences between the single and double gated 
structures. The paper is organized as follows: Section |l| 
contains a brief description of the theoretical formulation 
and the simulation procedure. Section [II presents the sim- 



ulation results along with an explanation of the underlying 
physics. Section IV summarizes this paper. 



II. Theory 

A schematic figure of a molecule coupled to gold contacts 
(source and drain) is shown in Fig. |l]a. As an example we 
use the Phenyl Dithiol (PDT) molecule which consists of 
a phenyl ring with thiol (-SH) end groups. A gate termi- 
nal modulates the conductance of the molecule. We use 
a simple model Hamiltonian H to describe the molecule 
(Fig. |l|b). The effect of the source and drain contacts is 
taken into account using self-energy functions Ei and E2 
Scattering processes may be described using another 
self-energy matrix Ep. However, in this paper we focus on 
coherent or ballistic transport (Ep = 0). The source and 
drain contacts are identified with their respective Fermi lev- 
els fii and iJ,2- Our simulation consists of iteratively solving 
a set of coupled equations (Fig. |l|c) - the Non-Equilibrium 
Green's Function (NEGF) formalism [0, 10 equations 
for the density matrix p and the Poisson's equation for the 
self-consistent potential Use- Given H, Use, ^ii E2, fii 
and 112 the NEGF formalism has clear prescriptions to ob- 
tain the density matrix p from which the electron density 
and the current may be calculated. Once the electron den- 
sity is calculated we solve the Poisson's equation to obtain 
the self-consistent potential Use- The solution procedure 
thus consists of two iterative steps: 

• Step 1: calculate p given Use using NEGF 

• Step 2: calculate Use given p using Poisson's equation 
The above two steps are repeated till neither Use nor p 
changes from iteration to iteration. It is worth noting that 
the self-consistent potential obtained by solving Poisson's 
equation (Eq. [ll] ) may be augmented by an appropriate 
exchange-correlation potential that accounts for many elec- 
tron effects using schemes like Hartree-Fock (HE) or Den- 
sity Functional Theory (DFT) [Q. In this paper we do 
not consider the exchange-correlation effects. 

A. Step 1: To obtain p from Use 

The central issue in non-equilibrium statistical mechan- 
ics is to determine the density matrix p; once it is known, 
all quantities of interest (electron density, current etc.) can 

^ The authors are aware of one experimental claim (J.H. Schon et 
al., Nature 413, page 713, 2001) reporting superior molecular FET 
with a single gated geometry. This claim, however, has been strongly 
questioned (see article by R.F. Service in Science 298, page 31, 2002). 
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Fig. 2 

Comparison of the simple tt orbital based model with an ab 
initio model (density functional theory with atomic orbital 

BASIS set) FOR THE PDT MOLECULE. THE ENERGY LEVELS CAN BE 
DIVIDED IN TWO SETS: OCCUPIED LEVELS (ANALOGOUS TO THE 
VALENCE band) AND UNOCCUPIED LEVELS (ANALOGOUS TO THE 

CONDUCTION band). ThE ENERGY GAP (ANALOGOUS TO THE 
BANDGAp) is THE ENERGY DIFFERENCE BETWEEN THE HIGHEST 

Occupied Molecular Orbital (HOMO) level (analogous to 

THE TOP OF the VALENCE BAND) AND THE LOWEST UNOCCUPIED 

Molecular Orbital (LUMO) level (analogous to the bottom 
OF the conduction band). The simple tt model agrees very 
well with the ab initio calculation in both the energy gap 

(middle) and homo and LUMO wavefunctions (left and 
right). The energy levels obtained from the simple model 
were equally shifted in energy so as to make the homo 
level coincide with the ab initio homo level. 



be calculated. A good introductory discussion of the con- 
cept of density matrix may be found in p^ . To obtain 
the density matrix p from the self-consistent potential Use 
using the NEGF formalism, we need to know the Hamil- 
tonian H, the contact self-energy matrices Si_2 and the 
contact Fermi levels /ii,2- In this section we describe how 
we obtain these quantities, and then present a brief outline 
of the NEGF equations. 

Hamiltonian: We use a simple basis consisting of one Pz 
(or tt) orbital on each carbon and sulfur atom. It is well 
known that the PDT molecule has tt conjugation - a cloud 
of TT electrons above and below the plane of the molecule 
that dictate the transport properties of the molecule p9[ |. 
The on-site energies of our pz orbitals correspond to the 
energies of valence atomic Pz orbitals of sulfur and carbon 
(apart from a constant shift of all levels which is allowed as 
it does not affect the wavefunctions). The carbon-carbon 
interaction energy is 2.5 eV which is widely used to de- 
scribe carbon nanotubes ]2C| ]. The sulfur-carbon coupling 
of 1.5 eV is empirically determined to obtain a good fit 
to the ab initio energy levels obtained using the commer- 
cially available quantum chemistry software Gaussian '98 
[0 (Fig. I)- 

Our model is very similar to the well established Pz or- 
bital based Hiickel theory used by many quantum chemists. 
Although we use a simple model Hamiltonian to describe 
the molecule, we believe that the essential qualitative 
physics and chemistry of the molecule is captured. This 
is because both the molecular energy levels and the wave- 
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The discrete levels of an isolated molecule broaden into a 

CONTINUOUS density OF STATES (DOS) UPON COUPLING TO 

CONTACTS. The gold FCC (111) contacts are modeled using a 

BASIS OF one 'S' type ORBITAL ON EACH GOLD ATOM. THE LUMO 
WAVEFUNCTION IS LOCALIZED ON THE PHENYL RING (FiG. ^ AND 
GIVES A SHARP PEAK IN THE DOS. THE HOMO IS DELOCALIZED AND 
GIVES A COMPARATIVELY BROADENED PEAK IN THE DOS. THE 
EQUILIBRIUM FeRMI ENERGY Ef (~ —5.1 eV FOR BULK GOLD) LIES 
■lUST ABOVE THE HOMO LEVEL. 



functions closely resemble those calculated from a much 
more sophisticated ab initio theory (Fig. ||). 

Self- energy: Self-energy formally arises out of partition- 
ing the molecule-contact system into a molecular subsys- 
tem and a contact subsystem. The contact self-energy E is 
calculated knowing the contact surface Green's function g 
and the coupling between the molecule and contact r. For 
a molecule coupled to two contacts (source and drain) the 
molecular Green's function at an energy E is then writ- 
ten as ]l6t (/: identity matrix, H: molecular Hamiltonian, 
Use- self-consistent potential): 

G=[EI-H -Use -^1-^2]-^ (1) 
where the contact self-energy matrices are 

Sl^2 = Ti^29l,2A.2 (2) 

We model the gold FCC (111) contacts using one s- 
type orbital on each gold atom. The coupling matrix el- 
ement between neighboring s orbitals is taken equal to 
—4.3 eV - this gives correct surface density of states (DOS) 
of 0.07 /{eV - atom) for the gold (111) surface ||. The 
site energy for each s orbital is assumed to be —8.74 eV 
in order to get the correct gold Fermi level of ~ —5.1 eV. 
The surface Green's function g is calculated using a re- 
cursive technique explained in detail in ]2^ ]. The contact- 
molecule coupling T is determined by the geometry of the 
contact-molecule bond. It is believed that when a thiol- 
terminated molecule like Phenyl Dithiol is brought close to 
a gold substrate, the sulfur bonds with three gold atoms 
arranged in an equilateral triangle. For a good contact ex- 
tended Hiickel theory predicts a coupling matrix element of 
about 2 eV between the sulfur pz orbital and the three gold 
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s orbitals. However to simulate the bad contacts typically 
observed in experiments |^ we reduce the coupling by 
a factor of five (this factor is also treated as a parameter, 
and our results do not change qualitatively for a range of 
values of this parameter). 

Unlike the Hamiltonian, the self-energy matrices are non- 
Hermitian. The anti-Hermitian part of the self-energy, also 
known as the broadening function: 



ri^2 = «(Si^2 - s| 2) 



(3) 



is related to the lifetime of an electron in a molecular eigen- 
state. Thus upon coupling to contacts, the molecular den- 
sity of states (Fig. ^ looks like a set of broadened peaks. 

Where is the Fermi energy?: When a molecule is cou- 
pled to contacts there is some charge transfer between 
the molecule and the contacts, and the contact-molecule- 
contact system attains equilibrium with one Fermi level 
Ef. A good question to ask is where Ef lies relative to the 
molecular energy levels. The answer is not clear yet, the 
position of Ef seems to depend on what contact model one 
uses. A jellium model 1^ for the contacts predicts that Ef 
is closer to the LUMO level for PDT whereas an extended 
Hiickel theory based model [E5| predicts that Ef is closer 
to the HOMO level (see Fig. g and the related caption for 
a description of HOMO and LUMO levels). Our ab ini- 
tio model [0 seems to suggest that for gold contacts, Ef 
(~ —5.1 eV) lies a few hundred millivolts above the PDT 
HOMO. In this paper we will use Ef = —5.1 eV and set the 
molecular HOMO level (obtained from the tt model) equal 

to the ab initio HOMO level ( 5.4 eV) (see Figs. § 1). 

Once the location of the equilibrium Fermi energy Ef is 
known, we can obtain the source and drain Fermi levels fix 
and 112 under non-equilibrium conditions (non-zero Vds)- 
Hi = Ef and 112 = Ef ~ qVos- 

NEGF equations: Given H, Si. 2, contact Fermi energies 
/ii^2 and the self-consistent potential Use, NEGF has clear 
prescriptions p6[ | to obtain the density matrix p. The den- 
sity matrix can be expressed as an energy integral over the 
correlation function —iG^{E), which can be viewed as an 
energy-resolved density matrix: 



dE[-^iG<{E)/2Ti] 



(4) 



The correlation function is obtained from the Green's func- 
tion G (eq. 1^) and the broadening functions Fi 2 (eq. H): 



iG< ^G{hTx + hT2)G^ 



(5) 



where fi^2{E) are the Fermi functions with electrochemical 
potentials 111^2'- 



ha{E) = 1 + exp 



E - Hi,2 



(6) 



The density matrix so obtained can be used to calculate 
the electron density n(r) in real space using the eigenvec- 
tors of the Hamiltonian ^'^(r) expressed in real space: 



iW = ^*«(r)*;(f)p„0 

a, 13 



(7) 



The total number of electrons N may be obtained from the 
density matrix as: 

N = trace(p) (8) 
The density matrix may also be used to obtain the termi- 



nal current |16|. For coherent transport, we can simplify 
the calculation of the current by using the transmission 
formalism where the transmission function [p^ : 

T{E) = trace [riGr2G^] (9) 

is used to calculate the terminal current 

/oo 
dET{E) {h{E)-f2{E)) (10) 
-00 

B. Step 2: To obtain Use from p 

The Poisson's equation relates the real space potential 
distribution U{r) in a system to the charge density n{r). 
We assume a nominal charge density no(r) obtained by 
solving the NEGF equations with U{r) — (at Vgs = 
Vds = 0). The Poisson's equation is then solved for the 
change in the charge density in — uq) from the nominal 
value rto : 



(eW(f)) =-q^{n{r)-no{r)) 



(11) 



The Poisson (or Hartree) potential U may be augmented 
by an appropriate exchange-correlation potential Uxc- In 
this paper, we do not take into account the exchange- 
correlation effects {Uxc — 0). We have two schemes to solve 
the Poisson's equation: (1) simple Capacitance Model and 
(2) rigorous numerical solution over a 2D grid in real space. 

Capacitance Model: We use a simplified picture of the 
molecule as a quantum dot with some nominal total charge 
A'o (at Vgs = Vds = 0) and some average potential U aris- 
ing because of the change N — Nq in this nominal charge 
due to the applied bias. Thus U, Nq and N are numbers 
and not matrices. The total charge N can be calculated 
from the NEGF density matrix using Eq. ^. U is the solu- 
tion to the Poisson's equation, and may be written as the 
sum of two terms: (1) A Laplace (or homogeneous) solution 
Ul with zero charge on the molecule but with applied bias 
and (2) A particular (or inhomogeneous) solution Up with 
zero bias but with charge present on the molecule. Thus 
U — Ul + Up. Ul is easily written down in terms of the ca- 
pacitative couplings Cms, Cmd and Cmg of the molecule 
(Fig. 0) with the source, drain and gate respectively: 



Ul = Pi-qVas) + ^^-^{-qVds) 



(12) 



where 



^ The potential distribution corresponding to the nominal charge 
density (when no drain or gate bias is anplied) is included in the 
calculation of the molecular Hamiltonian 
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Fig. 4 

Equivalent circuit model used to obtain the homogeneous 
(oh. zero charge) solution to the Poisson's equation. The 
molecular potential is controlled by the gate if the 

CAPACITATIVE COUPLING CmG BETWEEN THE MOLECULE AND THE 
GATE IS MUCH GREATER THAN THE CAPACITATIVE COUPLINGS CmS 
AND CmD BETWEEN THE MOLECULE AND THE SOURCE AND DRAIN 
RESPECTIVELY (SEE EQS. [12l [13 AND RELATED DISCUSSION). 
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MG 



c 



MS 



MD 



c 



MG 



(13) 



is a parameter (0 < /? < 1) and is a measure of how good 
the gate control is. Gate control is ideal when Cmg is very 
large as compared to Cms and Ca/d ■ In this case, (3=1 
and the Laplace solution Ul = —qVos is essentially tied to 
the gate. An estimate of gate control may be obtained from 
the numerical grid solution explained below by plotting [3 
as a function of gate oxide thickness (Fig. ||) . 

The particular solution Up may be written in terms of a 
charging energy Uo as: 



Up = Uo{N - No) 



(14) 



The charging energy is treated as a parameter, and may 
be estimated as follows. The capacitance of a sphere of 
radius R is iircR. If we distribute a charge of one electron 
on this sphere, the potential of the sphere is q/ineR. For 
R = 1 nm the value of this potential is about 1.4 eV. Thus 
we use a charging energy Uq 1 eV. Uq is the charging en- 
ergy per electron per molecule and may also be estimated 
from the numerical grid solution by finding the average po- 
tential in the region occupied by the molecule and carrying 
one electronic charge distributed equally. This numerical 
procedure also yields Uq ^ I eV and is used to estimate 
the charging energy while comparing the capacitance model 
with the numerical grid solution (see Fig. y and the related 
caption) . 

With the simple capacitance model just described, the 
Poisson's solution U is just a number. The self-consistent 

We have assumed that Cms = Cmd iu eq. which is reasonable 
because the center of the molecule is equidistant from the source and 
drain contacts in our model (see Fig. hi). In general, if the source 
(drain) is closer to the molecule, then Otis {Cmd) will be bigger Tl. 
With Cms = Cmd, the molecular Laplace poteiLtial is Vds/'^ in_tne 
absence of a gate (/3 = 0), as is evident from eq. [13 (also see Fig. te:;,d 
and the related caption). 



potential that adds to the Pz Hamiltonian (sec Eq. ^ is 
then calculated as Use = UI, where / is the identity ma- 
trix of the same size as that of the Hamiltonian. 

Numerical solution: We use a 2D real space grid to solve 
the discretized Poisson's equation for the geometry shown 
in Fig. 0a. The applied gate, source and drain voltages 
provide the boundary conditions. We use a dielectric con- 
stant of 3.9 for silicon dioxide and 2 for the self-assembled 
monolayer (SAM) |6|. 

The correct procedure to obtain the real space charge 
density n{r) (see Eq. ^l|) from the pz orbital space density 
matrix p is to make use of Eq. |. However, we simplify 
the calculation of n{r) by observing that a carbon or sulfur 
Pz orbital has a spread of about five to six Bohr radii (1 
Bohr radius ub — 0.529 A). So for each atomic site a we 
distribute a charge equal to paa equally in a cube with side 
~ lOap centered at site a. 

The solution to Poisson's equation yields the real space 
potential distribution. However, the self-consistent poten- 
tial UsG that needs to be added to the Hamiltonian (Eq. ^ 
is in the Pz orbital space. We assume that Use is a diag- 
onal matrix with each diagonal entry as the value of the 
real space solution U at the appropriate atomic position. 
For example, the diagonal entry in Use corresponding to 
the sulfur based Pz orbital would be equal to U{'rs) where 
fs is the position vector of the sulfur atom. 

111. Results 

The self-consistent procedure (Fig. ^) is done with the 
two types of Poisson solutions discussed above. The simple 
capacitance model is fast while the 2D numerical solution 
is slow but more accurate. The capacitance model has two 
parameters, namely /3 which is a measure of the gate con- 
trol, and Uq which is the charging energy. These parame- 
ters can be extracted using the 2D numerical solution. We 
will first present results with the capacitance model by as- 
suming ideal gate control, or f3 ~ 1. This ideal case is useful 
to explain the current saturation mechanism. We will then 
compare the results obtained from the capacitance model 
with those obtained from the numerical solution, and show 
that the two match reasonably well. 

A. Ideal gate control, on state 

Fig. H shows the molecular IV characteristic obtained by 
self-consistently solving the coupled NEGF - capacitance 
model Poisson's equations. We contrast the IV for ideal 
gate control {13 = 1, Fig. |^,b) with that for no gate control 
{(3 = 0, Fig. |c,d). For each case, we have shown the IV 
for positive as well as negative drain voltage. We observe 
the following: 

• With ideal gate control the IV is asymmetric with respect 
to the drain bias. For positive drain bias, we see very little 
gate modulation of the current. For negative drain bias we 
see current saturation and good gate modulation - the IV 
looks like that of a MOSFET. 

• With no gate control the IV is symmetric with respect 
to the drain bias. There is no gate modulation. 
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Fig. 5 

Three terminal molecular drain current vs. drain to source 

VOLTAGE CHARACTERISTIC WITH (a) IdEAL GATE CONTROL (/3 = 1), 
NEGATIVE DRAIN BIAS, (b) IdEAL GATE CONTROL, POSITIVE DRAIN 
BIAS, (C) No GATE CONTROL (/3 = 0), NEGATIVE DRAIN BIAS AND (d) 

no gate control, positive drain bias. with ideal gate 
control the iv is asymmetric with respect to drain bias. 
Good saturation and gate modulation is seen for negative 
drain bias but not for positive drain bias. with no gate 
control the iv is symmetric with respect to drain bias. for 
an explanation of the underlying mechanism for each of 
these iv curves, see fig. |6| 



These features of the IV characteristic may be under- 
stood as follows (Fig. Let us first consider the ideal 
gate case. Since the gate is held at a fixed potential with 
respect to the source, the molecular DOS does not shift 
relative to the source Fermi level as the drain bias is 
changed |^ . For negative drain bias (Fig. |^a), the drain 
Fermi level fj,2 moves up (towards the LUMO) with respect 
to the molecular DOS. Since the drain current depends on 
the DOS lying between the source and drain Fermi levels, 
the current saturates for increasing negative drain bias be- 
cause the tail of the DOS dies out as the drain Fermi level 
moves towards the LUMO. If the gate bias is now made 
more negative, the molecular levels shift up relative to the 
source Fermi level, thereby bringing in more DOS in the 
energy range between /ii and /i2 (referred to as the /Zi-/i2 
window from now on) , and the current increases. Thus we 
get current saturation and gate modulation. 

For positive drain bias (Fig. |b), fi2 moves down (to- 
wards the HOMO) with respect to the molecular DOS. The 
current increases with positive drain bias because more and 
more DOS is coming inside the fii-fj.2 window as 112 moves 
towards the HOMO peak. Once /i2 crosses the HOMO 
peak, the current levels off. This is the resonant tunneling 
mechanism. If the gate bias is now made more negative, 
no appreciable change is made in the DOS inside the /ii-/i2 




(c) 
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* This is true as long as the charging energy Uq 
typically the case. For high charging energies the 



Up can dominate the Laplace solution f/^ (see eqs 
discussion), thereby reducing gate control. 
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Fig. 6 

Gate induced current saturation mechanism: Assuming that 
the gate is very close to the molecule (ideal gate control, 
/3 = 1), the gate holds the molecular dos fixed relative to 

THE SOURCE FeRMI LEVEL fli BECAUSE THE GATE IS HELD AT A 

FIXED POTENTIAL WITH RESPECT TO THE SOURCE. WHEN A 
NEGATIVE DRAIN BIAS IS APPLIED (tOP LBFT), THE DRAIN FeRMI 
LEVEL ^2 MOVES UP RELATIVE TO THE MOLECULAR DOS. SiNCE THE 
DOS DIES OUT IN THE GAP, FOR SUFFICIENTLY HIGH DRAIN BIAS, NO 

MORE DOS COMES IN THE /il-/i2 WINDOW AND THE CURRENT 
SATURATES. WHEN A POSITIVE DRAIN BIAS IS APPLIED (tOP RIGHT) , 
IJ.2 MOVES DOWN RELATIVE TO DOS AND EVENTUALLY CROSSES THE 

HOMO. The IV is thus asymmetric. If the gate is far away 
(no gate control, 13 = 0), the DOS LIES roughly halfway 

between THE SOURCE AND DRAIN FeRMI LEVELS. In THIS CASE, FOR 

negative drain bias (bottom left), fj.1 CROSSES HOMO WHILE 
FOR POSITIVE DRAIN BIAS (BOTTOM RIGHT) fj,2 CROSSES HOMO. NO 
GATE MODULATION IS SEEN AS EXPECTED, AND THE CURRENT IS 
SYMMETRIC WITH RESPECT TO DRAIN BIAS. 



window, and the maximum current remains almost inde- 
pendent of the gate bias. 

Now let us contrast this with the case where no gate is 
present. In this case, due to the applied drain bias Vds, the 
molecular DOS floats up by roughly —qVosf^ with respect 
to the source Fermi level. For either negative (Fig. ^) or 
positive (Fig. |^d) drain bias, the current flow mechanism 
is resonant tunneling. Since the equilibrium Fermi energy 
lies closer to the HOMO, for negative drain bias crosses 
HOMO while for positive drain bias fi2 crosses HOMO , 
p7| . No gate modulation is seen as expected, and the IV 
is symmetric with respect to Vds- 

B. Ideal gate control - off state 

Fig. 1^ shows the log scale drain current as a function of 
gate bias at high drain bias. We note that despite assuming 
ideal gate control, the subthreshold slope of this molecular 
FET is about 300 mV/ decade which is much worse than 
the ideal room temperature ksT/q = 60 mV/ decade that 
a good MOSFET can come close to achieving. It is also 
worth noting here that our simulation is done at low tcm- 
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Fig. 7 

Subthreshold IV characteristic assuming ideal gate 
CONTROL. The temperature independent SUBTHRESHOLD SLOPE IS 
~ 300 mV/ decade EVEN WITH AN IDEAL GATE. This is because 

THE DOS IN THE HOMO-LUMO GAP DIES OUT VERY SLOWLY AS A 
FUNCTION OF ENERGY. THIS SLOW FALL OF THE DOS MAY BE 
ATTRIBUTED TO THE GOLD METAL-INDUCED GAP STATES. THUS A 
MOLECULAR FET WITH A RIGID MOLECULE ACTING AS THE CHANNEL 
IS A VERY POOR SWITCH. 



perature - the subthreshold slope of the molecular FET is 
temperature independent and only depends on the molecu- 
lar DOS as explained below. 

The poor subthreshold slope may be understood as fol- 
lows. As the gate voltage is made more positive, the molec- 
ular DOS shifts down with respect to the //i-/i2 window. 
The HOMO peak thus moves away from the ^i-/.t2 win- 
dow, and fewer states are available to carry the current. 
The rate at which the current decreases with increasing 
positive gate bias thus depends on the rate at which the 
tail of the DOS in the HOMO-LUMO gap dies away with 
increasing energy (Fig. ||). Typically we find that the tail 
of the DOS dies away at the rate of several hundred milli 
electron-volts of energy per decade, and this slow fall in the 
DOS determines the subthreshold slope. The slow fall in 
the molecular DOS may be attributed to the metal-induced 
gap (MIG) states - the gold source and drain contacts have 
a sizeable DOS near the Fermi energy, and are separated 
only by a few angstroms ^ Since the molecule is assumed 
to be rigid, the molecular DOS has no temperature de- 
pendence and hence neither does the subthreshold slope. 
Thus the molecular FET with a rigid molecule acting as 
the channel is a very poor switching device even with ideal 

^ For ballistic silicon MOSFETs, due to the band-limited nature of 
the doped silicon source/drain contacts, the MIG DOS is negligible. 
The subthreshold slope at a finite temperature is thus determined 
by the rate at which the difference in the source and drain Fermi 
function tails falls as a function of energy. This rate depends on 
the temperature, and the subthreshold slope is thus proportional to 
k«f/q (~ 60 mV at room temperature) for ballistic Si MOSFETs 
|14| . Preliminary results for a molecular FET with doped silicon 
source and drain contacts do show a subthreshold slope proportional 
to ksT/q; we are currently investigating this effect. 
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Fig. 8 

Estimate of gate control using 2D numerical Laplace's 

SOLUTION WITH A SINGLE GATED GEOMETRY (SOLID LINE) AND A 
DOUBLE GATED GEOMETRY (DASHED LINE). /3 (WHICH IS A MEASURE 
OF GATE CONTROL, SEE EQS. |i^, ^ AND RELATED DISCUSSION) IS 
PLOTTED AS A FUNCTION OF THE GATE OXIDE THICKNESS Tox- THE 
LENCJTH OF PDT MOLECULE (EQUAL TO THE CHANNEL LENGTH OF 
THE MOLECULAR FET) IS ABOUT 1 Jim. THUS IN ORDER TO GET 
GOOD GATE CONTROL (/3 > 0.8) THE GATE OXIDE THICKNESS HAS TO 
BE ABOUT ONE TENTH OF THE CHANNEL LENGTH, OR ABOUT 1 A! 
We have USED 3.9 and 2 as the dielectric constants FOR 
SILICON DIOXIDE AND THE SELF-ASSEMBLED MONOLAYER (SAM) 
RESPECTIVELY. FOR A DOUBLE GATED GEOMETRY, GOOD GATE 
CONTROL CAN BE OBTAINED WITH MORE REALISTIC OXIDE 
THICKNESSES (~ 7 A), AS EXPECTED SIMPLY BECAUSE TWO GATES 
CAN BETTER CONTROL THE CHANNEL THAN ONE. 



gate control. 

C. Estimate of Gate Control 
The 2D numerical Poisson's solution may be used to es- 



timate the gate control as follows, 
that 



From Eq. h2 we see 



f3 



1 dUL 



q dVc 



GS 



(15) 



Thus (3 may be estimated from the numerical solution by 
slightly changing Vqs (keeping Vds constant) and calcu- 
lating how much the Laplace's solution changes over the 
region occupied by the molecule. A plot of (3 calculated 
using this method as a function of the gate oxide thickness 
Tox is shown in Fig. |^. 

Knowing that the channel length (length of the PDT 
molecule) is about 1 nm. It is evident from Fig. || that in 
order to get good gate control (/3 > 0.8) the gate oxide 
thickness (Tox) needs to be about one tenth of the channel 
length {Lch), or about 1 A\ Thus we need Lch/Tox ^ 10 
to get a good molecular FET. In a well-designed conven- 
tional bulk MOSFET, Lch/Tox 40 [|l3|. This difference 
between a molecular FET and a conventional FET may be 
understood by noting that the dielectric constant of the 
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Comparison of the numerical Poisson solution with the 

CAPACITANCE MODEL. We SEE REASONABLE AGREEMENT IN THE TWO 
SOLUTIONS DESPITE THE SIMPLIFICATIONS MADE IN THE CAPACITANCE 
MODEL, (a) a REALISTIC CASE WITH tox = 1-5 nm WHICH YIELDS 

/3 = 0.28 AND Uo = 1.9 eV. No current saturation is seen, but 

SOME GATE MODULATION IS PRESENT, (b) An "ABSURD" CASE WITH 
tox = 1 A WHICH YIELDS /3 = 0.82 AND Uq = 1 eV (Uq IS LESS FOR 

THIS CASE BECAUSE THE GATE IS CLOSER TO THE MOLECULE; 
SCREENING EFFECT OF THE GATE ELECTRODE IS LARGER,). THE IV 
FOR THIS CASE LOOKS LIKE THAT FOR A MOSFET. SIMILAR IV MAY 
BE OBTAINED WITH tox = 1 nm, PROVIDED ONE USES A GATE 
INSULATOR WITH A DIELECTRIC CONSTANT ABOUT TEN TIMES THAT 
OF SILICON DIOXIDE. 



molecular environment 
that of silicon (==11.7) [ 



=2) is about 6 times smaller than 



Fig. H also shows as a function of Tox calculated using 
the 2D numerical Laplace's solution over a double gated 
molecular FET structure. In this case, we find that to get 
good gate control, we need Lch/Tox ~ 1-6. Thus a double 
gated structure is superior to a single gated one for a given 
Lch and Tox, as is also expected for conventional silicon 
MOSFETs. The reason for this is simply that two gates 
can better control the channel than one. 



D. Comparison: Capacitance model vs. Numerical Pois- 
son's solution 

Fig. H compares the IV characteristic obtained by solving 
the self-consistent NEGF-Poisson's equations with the nu- 
merical Poisson's solution and the capacitance model. The 
parameters f3 and Uq for the capacitance model were ex- 
tracted from the numerical solution. We see a reasonable 
agreement between the two solutions despite the simpli- 
fications made in the capacitance model (the capacitance 
model assumes a flat potential profile in the region occu- 
pied by the molecule, which may not be true, especially at 
high bias) . For tox = 1.5 nm (Fig. ^) there is very little 
gate modulation and no saturation as expected. In this 
case /3 = 0.28 (Fig. ||) and the IV resembles that shown in 
Fig. ||c more than the one in Fig. Also seen in Fig. || are 
the results for tox — 1 A. For this case /3 = 0.82 and we ob- 
serve reasonable saturation and gate control. For realistic 
oxide thicknesses, however, we expect to observe an IV like 
the one shown in Fig. We have also calculated the IV 
characteristics with a double gated geometry (not shown 
here), and as expected from Fig. ^, saturating IVs can be 
obtained for more realistic oxide thicknesses (~ 7 A). 

IV. Conclusion 

We have presented simulation results for a three terminal 
molecular device with a rigid molecule acting as the channel 
in a standard MOSFET configuration. The NEGF equa- 
tions for quantum transport are self-consistently solved 
with the Poisson's equation. We conclude the following: 

1. The current- voltage (IV) characteristics of molecular 
conductors are strongly influenced by the electrostatics, 
just like conventional semiconductors. With good gate con- 
trol, the IV characteristics will saturate for one polarity of 
the drain bias and increase monotonically if the polarity 
is reversed. By contrast two-terminal symmetric molecules 
typically show symmetric IV characteristics. 

2. The only advantage gained by using a molecular con- 
ductor for an FET channel is due to the reduced dielectric 
constant of the molecular environment. To get good gate 
control with a single gate the gate oxide thickness needs to 
be less than 10% of the channel (molecule) length, whereas 
in conventional MOSFETs the gate oxide thickness needs 
to be less than 3% of the channel length. With a double 
gated structure, the respective percentages are 60% and 
20%. 

3. Relatively poor subthreshold characteristics (a temper- 
ature independent subthreshold slope much larger than 
60 mV/ decade) are obtained even with good gate control, 
if metallic contacts (like gold) are used, because the metal- 
induced gap states in the channel preclude it from turning 
off abruptly. Preliminary results with a molecule coupled 
to doped silicon source and drain contacts, however, show a 
temperature dependent subthreshold slope {'^ kBT/q). We 
believe this is due to the band-limited nature of the silicon 
contacts, and we are currently investigating this effect. 

4. Overall this study suggests that superior saturation and 
subthreshold characteristics in a molecular FET can only 
arise from novel physics beyond that included in our model. 
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Further work on molecular transistors should try to exploit 
the additional degrees of freedom afforded by the "soft" (as 
opposed to rigid) nature of molecular conductors jl^ . 
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